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Abstract 

We describe results of Monte Carlo simulations on a model that seems to 
have the necessary ingredients to describe a disordered type-II superconductor 
in a magnetic field. We compute the free energy cost to twist the direction of 
the phase of the condensate and analyze the results by finite-size scaling. The 
results show convincingly that the model has different behavior as a function 
of dimension: in d = 4 the model clearly has a finite transition temperature; 
T c , while for d = 2 only there is only a transition at T = 0. 
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Since fluctuation effects play a much more important role in high temperature 
superconductors than in conventional superconducting materials, there has been a 
great deal of effort [0J to understand the behavior of type II superconductors in a 
magnetic field, including the effects of disorder, when one goes beyond the mean 
field picture of BCS or Ginzburg-Landau theories. One intriguing aspect which has 
emerged is the possibility of a vortex glass phase |2|, [3| in which the off diagonal long 
range order of the pair condensate has a phase which is random in space but frozen 
in time, much like the order parameter in a spin glass ||. This can arise because the 
Abrikosov flux lattice, which forms in pure samples, is destroyed by disorder in less 
than 4 dimensions |J beyond a certain length scale, i^ s . The phase of the condensate 
does not then form a regular periodic pattern on scales larger than ldi S , but, according 
to the vortex glass hypothesis, the system undergoes a transition into a spin glass- 
like state in which the phase is frozen in time. At the transition, the vortex glass 
correlation length, £ diverges. A number of experiments || have found evidence for 
such a transition in the I — V characteristics of Y-B-Cu-0 samples. Only if there is 
a vortex glass phase does the resistance really vanish M for H > H C1 . Otherwise, the 
resistance is, in principle, finite because clusters of vortices on scale £ can move by 
thermal activation over barriers, a process known as "flux creep" . These effects are 
observable || in high-T c compounds since they have much larger fluctuations than 
conventional materials. 

In this paper we descrive results of Monte Carlo simulations, analyzed by finite- 
size-scaling techniques, on a model that seems to have the necessary ingredients to 
describe the vortex glass state. The results show convincingly the different behavior 
occurs in different dimensions. In d = 4 there is clearly a transition at finite transition 
temperature T c with vortex glass order at lower temperature, while in d = 2 there is 
only a transition at T = 0. Analogous results for d = 3 have been presented before 
0, and indicate behavior close to what is expected at the lower critical dimension, 
though with probably a finite T c . Our results for d — 4 have been briefly described 
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in a conference proceeding JXO | . We feel that it is useful to include them here as a 
contrast to the new results for d = 2 to emphasise the power of the finite-size-scaling 
technique in elucidating whether or not a glass-like transition occurs. 

The model that we study, known as the "gauge glass", has the following Hamil- 
tonian: 

H = - ^2 COS (0i _ <t>j ~ A ij) ■ (!) 
<i,j> 

The phase, (pi, is defined on each site of a regular lattice, square for two dimensions, 
simple cubic for d = 3 and simple hypercubic for d = 4, with N = L d sites. Periodic 
boundary conditions are imposed. The sum is over all nearest neighbor pairs on the 
lattice. The effects of the magnetic field and disorder are represented by the quenched 
vector potentials, A^, which we take to be independent random variables with a 
uniform distribution between and 2ir. This model seems to be the simplest one with 
the correct ingredients of randomness, frustration and order parameter symmetry. It 
does, however, ignore screening, and therefore corresponds to an extreme type II limit 
in which k = A/£ — ► oo, where A is the penetration length. Since k > 1 in the high 
T c superconductors, this limit is not unreasonable. It is unclear, however, how much 
inclusion of screening via a fluctuating gauge field would modify the behavior of Eq. 

(0)- 

If the Aij are restricted to the values and tt, the model becomes the XY spin 



glass, for which the lower critical dimension is believed [11 1 to be 4. However, earlier 



work ]T3L |T^, ||, has shown that the gauge glass is in a different universality class 
from the XY spin glass, presumably because it does not have the the "reflection" 
symmetry, (pi — > —(pi Vz |13| . 

As discussed before || [T2], it is useful to consider the change in free energy 
AF when one imposes a twist G along one of the space directions, x say. More 
precisely, the periodic boundary conditions, (pi = (pi+ix are replaced by the twisted 
boundary conditions, (pi = (pi+Lx + @- By a simple redefinition of the phases (pi one 
can replace this situation by a system with periodic boundary conditions and an extra 
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contribution, Q/L, to the vector potential on bonds in the x-direction. 

By Monte Carlo methods one can calculate derivatives of the free energy w.r.t. 
6, so, for a single sample, we define a current, /, and a stiffness, Y, by 



r = ^ = ^{E( cosA ^)r-^E[( sinA ^ inA J )T- (sinA i ) T (sinA J ) T ]} , (3) 

where A, = 0« — — A iti+ x, F is the total free energy and i + x refers to the nearest 
neighbor site in the x-direction from i. Note that both / and Y are gauge invariant 
so they are still useful even if one includes fluctuating gauge fields. 

Above T c , AF, and hence both I and Y, go to zero rapidly with increasing system 
size because the system is insensitive to boundary conditions when L is much greater 
than the vortex glass correlation length £. If T c is finite, then, below T c , I and Y 
vary with L as L 6 where 9 (> 0), is an exponent describing the low temperature 
phase. In other words, I and Y increase with increasing L below T c , the opposite 
of what happens above T c . Precisely at T c , both I and Y are independent of size. 
Hence if T c is finite, I and Y should come together at T c and splay out again at lower 
temperatures. By contrast, if T c — 0, then, at T — 0, I and Y vary as L e but with 
9 < 0. Consequently, / and Y decrease with L even at T = 0. 

In a disordered system, it is necessary to perform an average over different real- 
izations of the disorder, which we indicate by [• • -] av . For the gauge glass, the average 
values of / and Y are both zero, i.e. 

[Y] av = [I] av = , (4) 

because the configuration in which the vector potentials in the x direction have been 
increased by Q/L has the same weight in the configurational average as the original 
choice of vector potentials. One is therefore interested in the root mean square fluc- 
tutation between samples. This means that many samples must be averaged over, 
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typically several thousand. If T c is finite, the finite size scaling form for the r.m.s. 
current, AI is therefore 



AI = [I 2 ] V 2 = I(L^(T - T c )) (T c >0) , (5) 

where v is the correlation length exponent. We shall concentrate on the r.m.s. current 
in what follows, rather than the stiffness, because sample to sample flucutations in 
the stiffness have an asymmetric distribution with a long tail, which makes it difficult 
to get good statistics [PJ. If T c — then AI decreases with size even at T = 0, i.e. 
AI ~ L 9 where 9 is negative and related to the exponent v giving the divergence of 
the correlation length as T — > by —9 = \jv. The finite size scaling form is then 

L 1/U AI = I(L 1/v T) (T c = 0) . (6) 



Tests to ensure equilibration were carried out as described elsewhere [|4 



We first show results for d = 4, [10]. Fig. 1 shows clearly that the data for 



the r.m.s. current for different sizes come together at T = T c ~ 0.95 and then 
splay out again on the low temperature side, just as expected at a finite temperature 
transition. This provides unambigous evidence that there is a vortex glass transition 
in four dimensions (and presumably also in higher dimensions) with vortex-glass order 
on the low-temperature side. A scaling plot corresponding to Eq. (||) is shown in Fig. 
2. From the fit we estimate 

T c = 0.96 ±0.01, i/ = 0.7 ±0.15 (d = 4) . (7) 

Next we discuss the case of d — 2. The results for AI are shown in Fig. 3. Notice 
that they are quite different from Fig. 1, since, even at the lowest temperature, AI 
decreases with increasing size. This is precisely what is expected at a zero temperature 
transition, and the scaling plot in Fig. 4 corresponding to Eq. (|6|) works very well. 
From the fit we estimate 

T c = 0, v = 2.2 ± 0.2 (d = 2) . (8) 
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The value for v agrees with earlier work |T2 |, in which a different finite-size-scaling 
technique was used. Recent experiments on very thin (16A) films of YBCO have 
provided striking confirmation that T c = for the vortex glass in two dimensions. 
Furthermore, it is found that non-linear current-voltage characteristics set in when 



the current density exceeds a value, J n \ which, according to scaling theory [12j varies 
with temperature as T l+U . According to conventional flux creep theory, there is no 
divergent length scale as T — > 0, which corresponds to setting v = 0. The experiments 



15| find 1+v = 3.0±0.3 in excellent agreement with the results presented here and in 



Ref. flzfl ) but in clear disagreement with the flux creep theory. We emphasize, then, 
that the vortex glass picture leads to measurable consequences even when T c = 0. 

To conclude, we have shown that a Monte Carlo calculation of the change in free 
energy due to a twist in the boundary conditions combined with finite size scaling 
is a very powerful tool for systems with XY-like symmetry. We have shown that it 
clearly distiguishes between the finite temperature transition in the d = 4 gauge glass 
and the zero temperature transition of the two-dimensional model. 
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Figure Captions 

Figure 1: The r.m.s. current, A/ = [I 2 ]^ 2 , for d = 4 determined by Monte Carlo 
simulations for different sizes and temperatures. The curves for different sizes 
are expected to come together at T c and, if there is order in the low temperature 
state, to splay out again at lower temperatures. The data does indeed follow 
this behavior. 

Figure 2: The same data as in Figure 1 but in a finite size scaling plot, with 
T c = 0.96 and v = 0.7. 

Figure 3: The r.m.s. current, AI = [I 2 ]^ 2 , for d = 2 determined by Monte Carlo 
simulations for different sizes and temperatures. The curves for different sizes 
do not come together, even at the lowest temperature. This behavior indicates 
a transition at T = 0. 

Figure 4: The same data as in Figure 4 but in a finite size scaling plot, with T c = 
and v = 2.2. 
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Figure 1 
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